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We apply the energy surface method to study a system of Na three-level atoms interacting with 
a one-mode radiation field in the A- and ^/-configurations. We obtain an estimation of the 
ground- state energy, the expectation value of the total number of excitations, and the separatrix of 
the model in the interaction parameter space, and compare the results with the exact solutions. We 
have first- and second-order phase transitions, except for the ^/-configuration which only presents 
second-order phase-transitions. 
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1. Introduction. The Tavis-Cummings Model [1], 
which describes the interaction of a collection of N two- 
level atoms with a quantized electromagnetic field in the 
dipolar and rotating- wave approximations (RWA), has 
an extensive use in quantum optics [2]. Recently this 
model has been physically realized using a QED cavity 
with Bose-Einstein condensates [3, 4]. Particularly inter- 
esting has been the investigation of the phase transitions 
of the system in the thermodynamic limit [5] , and at zero 
temperature [6, 7]. 

The physics of three-level systems interacting with one 
or two quantized modes of the electromagnetic field is 
very rich and many special dynamical situations have 
been studied. In particular, a formalism to describe one 
three-level atom interacting with a one- or a two-mode 
field has been discussed together with the atomic level 
occupation probabilities, coherence properties, photon 
probability distribution, fluctuations, and the evolution 
of squeezing in a series of works [8]. For one three-level 
atom interacting with a one- or a two-mode field, it was 
found that the phase distribution properties of the field 
reflect the collapses and revivals of the level occupation 
probabilities. However, for the two- mode case, there are 
exceptions and the collapses and revivals are decorrelated 
from the phase field [9]. 

A comprehensive review of the dynamical interac- 
tion of an atom with radiation in the framework of 
the Jaynes-Cummings type has been done by Yoo and 
Eberly [10]. However, by means of the thermodynamic 
Green function the phase transitions of the Dicke model, 
including all modes of the radiation field, for atoms con- 
fined in a cubic resonance wavelength, have been consid- 
ered for a finite temperature in [11]. Recently, there has 
been a semi-classical treatment of the phase transitions 
in a system of three-level atoms in the A-configuration 
interacting with a two-mode quantized electromagnetic 
field without the RWA, using however the generalized 
Holstein-Primakoff transformation to find the separatrix 
of the system in the thermodynamic limit [12]. 

The main contribution of this manuscript is the es- 
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FIG. 1. Schematic representation of the three atomic config- 
urations, huji denotes the energy of the i-th level, and /j^ij 
the dipolar coupling between levels i and j. We use 7i = 1 
throughout. 



tablishment of the separatrix of the three-level system 
interacting with a one-mode electromagnetic field for the 
three atomic configurations, S, A, and V (cf. Fig. 1). 
Each separatrix determines in control parameter space 
(dipolar strenghts) where the quantum phase transitions 
take place, and which are of first- and second-order de- 
pending on the values of the dipolar interactions. The 
presence of the quantum phase transitions can be clearly 
seen in the calculation of the ground state energies and in 
the expectation value of the total number of excitations, 
as functions of the dipolar couplings /ii2, /ii3, and /i23- 
The agreement with the corresponding exact quantum 
calculations is remarkable, in spite of considering a small 
number Na of atoms {Na = 2 and Na = 10). Further- 
more, by taking /i23 = in the E- and A- configurations, 
and /ii3 = in the V case, the corresponding Hamilto- 
nian systems describe two-level problems in the manner 
of Tavis-Cummings, giving consistency to our results. In 
these cases it is straightforward to check that the sep- 
aratrix coincides with the one established by Hepp and 
Lieb [5]. 

The organization of the paper is as follows. In Section 
2 we establish the model and present the basic formula- 
tion of the problem. In Section 3 an analytic expression 
for the energy surface of the system is obtained, and use 
it to obtain an estimation of the ground state of the sys- 
tem in the three basic configurations. In Section 4 we 
solve the Hamiltonian for a finite number of particles 
and compare the results with those previously obtained. 
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Finally, a summary of the general results and conclusions 
is outlined in Section 5. 

2. Model Hamiltonian. We consider a quantum system 
of Na three-level atoms, each atom being able to occupy 
one of three levels characterized by energies cji, UJ2 and 
cja, interacting dipolarly with a one-mode field of fre- 
quency ^1, and assume the rotating- wave approximation. 
We consider uji < uj2 < ^3- The Hamiltonian describing 
this system can be written as [10] 



k=i 



Ml3 



M23 



(a iai + at A,s) - ^ ^2 + at ^3) (1) 



where the first two terms on the r.h.s. are related to the 
free atomic and field parts. The a, at, denote the annihi- 
lation and creation operators of the radiation field, and 
the operators Aij are the generators of a U{3) algebra. 
The operators a, at, satisfy the Heisenberg-Weyl algebra 
commutators 



[a, at] 



a, a 







while the atomic operators Aij (weight for z = j, lowering 
for i > j and raising with i < j) fulfill the [/(3)-algebra 
commutation relations 



Aij^ Aj^i 



SjkAii - SiiA 



i^kj 



The usual three-level atomic arrangements are obtained 
from (1) imposing /iis = for the S-configuration, 
M12 = for the A-configuration, or /i23 = for the V- 
configuration. In each case there are two constants of 
motion: the total number of atoms 



Na = All + A22 + As3 , 

and the total number of excitations 

Ms = ata + i22 + 2i33 



(2) 



Ma a^a + A33 

^22 ■ 



My ata 



■A 



33 
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3. Energy surface. In order to obtain an energy sur- 
face, we use as a trial state the direct product of coherent 
states in each subspace. The use of coherent states as 
trial states lets us determine in analytic form the expec- 
tation values of matter and field observables. Thus we 
use the tensor product of Heisenberg-Weyl HW{1) co- 
herent states for the radiation part, \a} = |0), and 
U{3) coherent states for the atomic part. The unnormal- 
ized U{3) coherent state can be constructed by taking 
the exponential of the lowering generators acting on the 
highest weight states of /7(3), 



|[^i,^2,/^3]7i^72,73} = e 



73^21 ^72^31 ^71^32 1 



[hi,h2,hs]) , 




where we denote by 
I [^1, ^2, ^^3]) 



the highest weight state of the Gelfand-Tsetlin basis of 
the irreducible representation [/ii, ^3] of U{3) [13]. 
The action of the raising operators Aij on the highest 
weight state vanishes: Aij \ [hi, ^3]) = 0. 

In this contribution we only consider the completely 
symmetric representation [Na, 0, 0] of U{3), where Na 
denotes the number of atoms. Therefore, the trial state 
is the tensor product 



|O)0 e^'^^' e^'^''\[Na,0,0]) . 



The parameter 71 does not appear in the coherent state 
because A32 \ [Na, 0, 0]) = 0. 

The expectation value with respect to this state of the 
model Hamiltonian (1) gives the energy surface 



^(^, ^2, ^3, ^1, ^2, ^3) = ^g'^ 

2^aQ 



■N 



Ul + CJ2 ^3 + ^3 Q2 



1 + ^2 + d 
M12 ^3 COS ^3 + Ml3 Q2 COS ^2 + /i23 Q2 QS COS l^i 



1 + ^2 + d 

where we have used the identifications = Qk exp(z cpk), 
with k = 2, 3, and a = gex.p{i(p). Additionally, the 
angles = - (^2 + ^3, ^2 = ^ - ^2, and ^3 = - (^3 
were defined. 

The minima for the energy surface are obtained when 
'i^i = 0, TT for z = 1, 2, 3, and 

/il2 COS'i^3 > 0, /ii3 COS'i^2 > 0, j^23 COS'di > . 

To have an intensive quantity we divide the energy sur- 
face 1-L{q, Q2-, Q?>-,^i-,^2-,^?>) by Na, and use the conditions 
of a minimum for the angles ^i to obtain 



^(^, ^2, ^3) = Q 
-2-Q 



.2 , UJl + CJ2 ^3 + ^3 Q2 
M12I ^3 + |M13| Q2 + 1^23 1 Q2 QS 



1 + ^2 + ^3 



(4) 



where p = p/^/N^, and uok and pki are now measured in 
units oi Vt. In what follows we will set 1] = 1. 

Using the Ritz variational principle, one finds the best 
variational approximation to the ground-state energy of 
the system and its corresponding eigenstate for the en- 
ergy surface (4). The values of Qc are always given in 
terms of Q2c and ^3c. These last two, independent of 
Na, are found numerically. For the configuration S the 
estimation of the ground energy is plotted in Fig. 2 (a) 
as a function of pi2 and /i23- In Fig. 2 (b) we show also 
the corresponding value for the constant of motion M^. 
In both figures we used the resonant case CJ32 =^0^21 = 1, 
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FIG. 2. (color online) Numerical results for the lowest value 
of the energy surface for the S-configuration (a) and its cor- 
responding estimation for the constant of motion Ms (b) , 
as functions of /ii2 and //23- We consider the resonant case 
UJ32 — 0J21 — 1, with cji = 0. 




FIG. 3. (color online) Numerical results for the lowest value 
of the energy surface for the A-configuration (a) and its cor- 
responding estimation for the constant of motion Ma (b) , as 
functions of /iis and /i23. We consider the non-resonant case 
cjai — 1.3, CJ32 — 0.8, and uji — 0. 



and cji = 0, with the notation ujij = uji — ujj. The white 
lines indicate the separatrix of the configuration. 

For the configuration A the result for the semiclassical 
ground energy is plotted in Fig. 3 (a) as a function of 
/ii3 and /i23- In Fig- 3 (b) we show the corresponding 
value of the constant of motion M\. In both figures we 
used ujsi = 1-3, 6^32 = 0.8, and uji = 0. Note that, in this 
case, we chose for illustrative purposes to work away from 
resonance. Again the white lines denote the separatrix 
of the system. 

Finally, considering the V-configuration we obtain the 
semiclassical ground energy shown in Fig. 4 (a) as a func- 
tion of /ii2 and /ii3, and in Fig. 4 (b) the corresponding 
value of the constant of motion My is displayed. In both 
figures we used again the resonant case ujsi = UJ21 = 1, 
with uji = 0. The white line is the corresponding sepa- 
ratrix for the configuration. 

The separatrix of the system is obtained by analizing 
the stability and equilibrium properties of the energy sur- 
face by means of the catastrophe formalism. Thus, the 
locus of the points where the thermodynamic phase tran- 
sition M = ^ M 7^ occurs, is given in each case by: 
for the S-configuration 

1^12 + [|M23| - ^A^f B (IM23I - V^) = ^21 ; (5) 
for the A-configuration 

M?3 + [Im23| - (IM23I - = ^31 ; (6) 




FIG. 4. (color online) Numerical results for the lowest value 
of the energy surface for the ^-configuration (a) and its cor- 
responding estimation for the constant of motion My (b) , as 
functions of /ii2 and /ii3. We consider U31 — U21 — 1 with 
ui 0. 

and for the V-configuration 

^ + ^ = 1 . (7) 

CJ2I ^31 

For the S-configuration, the phase transition across 
M12 = \/^2i in the separatrix is of second-order, and the 
one that takes place along the segment of the circum- 
ference is of first-order. For the A-configuration some- 
thing similar happens: across /ii3 = y^cJsi in the sepa- 
ratrix the phase transition is of second-order, and along 
the segment of circumference in the separatrix is of first- 
order. For the V-configuration all transitions are always 
of second-order. 

4. Quantum case. In the quantum case we must con- 
struct the proper basis for each configuration taking into 
account (3). Thus, the basis is characterized by A/'^, the 
number of atoms, which besides defines the U (3) highest 
weight state, and M, the total excitation number for the 
configuration considered. 

The value of the matrix elements of the U (3) genera- 
tors in the Hamiltonian (1) in the Gelfand-Tsetlin basis 
are given in [14]. Using the matrix elements, we calcu- 
late the matrix Hamiltonian and through a diagonaliza- 
tion procedure obtain the ground state energies, and the 
expectation value of the total number of excitations. In 
some cases for all the configurations one can get analytic 
expressions for the lowest energy states for a few number 
of atoms. For example, for the S-configuration and one 
atom we obtain 

which coincides with the results given in [10] after an 
identification of the parameters used here. 

In Figs. (5), (6), (7), we show the ground-state energy 
and the value of the constant of motion M for the 
A-, and V-configurations as functions of the interaction 
parameters for Na = 2 atoms, and for the same param- 
eter values as in the semiclassical case, for comparison 
purposes. We can observe that the resemblance between 
the semiclassical and quantum results is excellent. In the 
quantum case we also take notice of phase crossovers in 
the ground-state energy which occur every time the value 
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FIG. 5. (color online) Numerical results for the ground-state 
energy of the model Hamiltonian for the S-configuration (a) 
and the corresponding value of the constant of motion Ms (b) , 
as functions of /J12 and /i23. We consider the case Na = 2, 
UJ32 — 0021 — 1 and cji = 0. The white line corresponds to the 
semiclassical separatrix. 




FIG. 6. (color online) Numerical results for the ground-state 
energy of the model Hamiltonian for the A-configuration (a) 
and the corresponding value of the constant of motion Ma (b) , 
as functions of /xis and /i23. We consider the case Na — 2, 
— Q = 0.3, UJ32 — ^ = —0.2 and uji — 0. The white line 
corresponds to the semiclassical separatrix. 

of the constant of motion changes [15]. This may be seen 
as different shades in the energy surface, and is better 
seen in the plot for M where there is a jump when a 
phase crossover takes place. For the S-configuration we 
see in Fig. (5) that there are phase transitions from the 
region Ms = directly to regions Ms = 1, 2, . . . , 5, for 
Na = 2 atoms. When the number of atoms increases, 
greater values of Ms are pulled towards the region with 
Ms = 0, in such a way that there will be phase tran- 
sitions from this value to greater Ms's. For the other 
configurations, one has a similar effect when the number 
of atoms increase. 

When the number of atoms Na increases, all the curves 
where the constant of motion changes tend to the classi- 




FIG. 7. (color online) Numerical results for the ground-state 
energy of the model Hamiltonian for the ^/-configuration (a) 
and the corresponding value of the constant of motion My 
(b), as functions of /J12 and /iis. We consider the case Na = 2, 
(^31 = = 1 and cji = 0. The white line corresponds to the 
semiclassical separatrix. 




FIG. 8. (color online) Comparison of the equipotential curves 
of the ground-state energy of the model Hamiltonian in the 
S-configuration for Na = 2, 10, lower continuous and middle 
broken lines respectively, with CJ32 = CJ21 = 1 and cji = 0. 
The upper line corresponds to the semiclassical separatrix, to 
which all tend as A^a ^ 00. 

cal separatrix. This behavior is shown in Fig. (8), where 
we plot the equipotential curves for the ground-state en- 
ergy of the model in the S-configuration with Na = 2, 10 
atoms for Ms = 0, and compare with the semiclassical 
separatrix. 

5. Conclusions. We studied the three main configura- 
tions of a system of Na three-level atoms interacting with 
a one-mode radiation field, which is the simplest general- 
ization of the Tavis-Cummings model. Using HW{1)- 
and 5'[/(3)-coherent states, we established the energy 
surface for the system and obtained an approximate ex- 
pression for the ground-state of the system, as well as the 
separatrix in the parameters of the model which defines 
the locus of the quantum phase transitions. Although 
the phase transitions in all the configurations are simi- 
lar, in the thermodynamic limit there are only two zones, 
and the precise form for the V-configuration is qualita- 
tively different to the S- and A-configurations. For the 
S- and A-configurations the phase transition from one 
region to the other can be of first- or second-order, de- 
pending on the zone where the separatrix is crossed; for 
the V-configuration all phase transitions are of second- 
order. The comparison with the exact ground-state of the 
model shows that our approximation is excellent, at least 
for the energy and the expectation value of the constant 
of motion of each configuration. We were also able to 
verify numerically that when the number of atoms goes 
to infinity, the multitude of quantum phase crossovers 
tend to the thermodynamic limit. 
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